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ABSTRACT 



X 



Aims. We present for the first time a coherent model of the polarized Galactic synchrotron and 
thermal dust emissions which are the main diffuse foreground for the measurement of the po- 
larized power spectra of the CMB fluctuations with the Planck satellite mission. 
Methods. We produce 3D models of the Galactic magnetic field including regular and turbulent 
components, and of the distribution of matter in the Galaxy, relativistic electrons and dust 
grains. By integrating along the line of sight we construct maps of the polarized Galactic syn- 
chrotron and thermal dust emission for each of these models and compare them to currently 
available data. We consider the 408 MHz all-sky continuum survey, the 23 GHz band of the 
Wilkinson Microwave Anisotropy Probe and the 353 GHz Archeops data. 

Results. The best-fit parameters obtained are consistent with previous estimates in the liter- 
ature based only on synchrotron emission and pulsar rotation measurements. They allows us 
to reproduce the large scale structures observed on the data. Poorly understood local Galactic 
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structures and turbulence make difficult an accurate reconstruction of the observations in the 
Galactic plane. 

Conclusions. Finally, using the best-fit model we are able to estimate the expected polarized 
foreground contamination at the Planck frequency bands. For the CMB bands, 70, 100, 143 and 
217 GHz, at high Galactic latitudes although the CMB signal dominates in general, a significant 
foreground contribution is expected at large angular scales. In particular, this contribution will 
dominate the CMB signal for the B modes expected from realistic models of a background of 
primordial gravitational waves. 

Key words. ISM: general - Methods: data analysis - Cosmology: observations - Submillimeter 
1. Introduction 

The PLANCK satellite mission, currently in flight, should permit the more accurate 
measurements of the CMB anisotropics both in temperature and polarization. Planck, 
which observes the sky on a wide range of frequency bands from 30 to 857 GHz, 
has a combined sensitivity of ~ and an angular resolution from 33 to 5 

arcmin (jPLANCK b lucboo k 2004)) . Of particular interest is the measurement of the 
polarization B modes which implies the presence of tensor fluctuations from primordial 
gravitational waves generated during inflation. Planck should be able to measure 
the tensor-to-scalar ratio, r, down to 0.1 ([Bctoule ct al 2009, Efstat hiou et al 2009|l 
in the case of a nominal mission (2 full-sky surveys) and to 0.05 for the Extended 
Planck Mission: four full-sky surveys (Efstathiou fc Gratton 2009[) . The value of r sets 
the energy scale of the inflation (|Peiris et al 2003]) and then provides constraints on 
inflationnary models ([Baumann et al 2009[) . 

To achieve this high level of sensitiviy it is necessary to accurately estimate 
the temperature and polarization foregrounds mainly from diffuse Galactic emission 
components - synchrotron, thermal and rotational dust, and free-free - as well as from 
the Galactic and extra-Galactic point-like and compact sources. Indeed, at the Planck 
frequency bands these foreground components may dominate the CMB signal and 
therefore, they need to be either masked or subtracted prior to any CMB analysis. For 
this purpose, the Planck collaboration plans to use component separation techniques 
(see ILeach et al 20081 for a summary) in addition to the traditional masking of highly 
contaminated sky regions including identified point-like and compact sources. As these 
component separation techniques will be mainly based on Planck data only, one of the 
main issues will be to estimate the residual foreground contamination on the final CMB 
temperature and polarization maps. These residuals will translate into systematic biases 
and larger error bars on the estimation of the temperature and polarization power 
spectra of the CMB fluctuations (see lBetoule et al 20091 for a recent study). Thus, they 
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will impact the precission to which cosmological information can be retrieved from the 
Planck data. 



For polarization the main foreground contributions will come from the diffuse 
Galactic synchrotron and thermal dust emission. From the WMAP observations, 
Page et al 2007| have shown that the radio synchrotron emission from relativistic 
electrons is highly polarized, up to 70 %, between 23 and 94GHz. Furthermore, 
IBenoit et al 20041 IPonthieu et al 20051 have observed signicantly polarized thermal 
dust emission, up to 15 % at the 353 GHz Archeops channel. By contrast the diffuse 
free-free emision is not polarized and the anomalous microwave emission has been 
measured to be weakly polarized, 3^'g %, (|Battistelli et al 2006|l . Finally, at the 
Planck frequency bands the polarized contribution from compact and point sources is 
expected to be weak for radio (jNolta et al 2009|) and dust (jDesert et al 2008)) sources. 
The spatial and frequency distribution of both Galactic synchrotron and thermal dust 
polarized emission at the Planck frequencies are not well known and the only available 
information comes from micowave and submillimctrc observations. For synchrotron, 
Faraday rotation (jBurn et al 1966[) makes it very difficult to extrapolate the polarized 
observed radio emission (jWolleben et al 20061 IWolleben et al 20071 ICarretti et al 2009j) 
to the microwave domain. For thermal dust, polarized observations are not currently 
available in the infrared and the current optical data (JHeilcs 2000]) are too sparse 
( |Page et al 2007[ ) for a reliable extrapolation to lower frequencies. 



The diffuse Galactic synchrotron emission is produced by relativistic electrons 
spiraling around the Galactic magnetic field lines and its polarization is orthogonal 
both to the line-of-sight and to the field direction (Ribicki & Lightman 1979). Based on 



these statements, Page et al 2007 proposed to model the polarized synchrotron Galactic 
emission observed by the WMAP satellite using a 3D model of the Galaxy including the 
distribution of relativistic electrons and the Galactic magnetic field structure. Although 
this model allowed them to explain the observed polarization angle at the 23 GHz band 
where the synchrotron emission dominates, it was not used for the CMB analysis. Instead 
they constructed a template of the polarized synchrotron emission from the 23 GHz band 
and extrapolated it to higher frequencies. Independently. Ban et al 20041 IHan et al 2 006 
has also proposed a 3D model of the Galactic free electrons (|Cordes & Lazi o 2002]) 
and of the Galactic magnetic field including a regular and a turbulent component to 
explain the observed rotation measurements on known pulsars. Based on previous works 
ISun et al 20081 has performed a combined analysis of the polarized WMAP data and 
of the rotation measurements of pulsars using the publicly available HAMMURABI 
code (Waolkcn s et al 2009[) for computing the integrated emission along th line-of-sight. 
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This work has been extended by IJaffe et al 20091 for the study of the Galactic plane 
using a MCMC algorithm for the determination of the parameter of the models and 
by Jansso n~et al 200"9l for the full sky using a likelihood analysis for parameter estimation. 

Dust grains in the Interstellar Medium (ISM) are heated by stellar radiation and 
radiate in the form of thermal dust emission (jDesert et al 1998|) . They are considered to 
be oblate and to align with their longitudinal axis perpendicular to the magnetic field 
lines (jDavis fc Greenst cin 195lj) . When aligned they will end up rotating with their 
angular moment parallel to the magnetic field direction. As the thermal dust emission 
is more efficient along the long axis, this generates linear polarization orthogonal to the 
magnetic field direction and to the line-of-sight. The polarization fraction of the emission 
depends on the distribution of the size of the grains and is about a few percent at the 
millimeter wavelengths (jHildebrand et al 19991 IVaillancourt 2002 j) . IPonthieu et al 20051 
concluded that the polarized emission observed in the 353 GHz Archeops data was 
associated with the thermal dust emission and proposed a simple magnetic field pattern 
to explain the measured direction of polarization on the Galactic plane. |Page et al 2 007 
suggested that part of the observed polarized emission at the 94 GHz WMAP data was 
also due to thermal dust. They modeled it using the observed polarization direction of 
stellar light (Heilcs 2000) which from the above statement must be perpendicular to 
that of thermal dust. 

With the prospect of the data analysis of the Planck satellite mission in mind, we 
present here consistent physical models of the synchrotron and thermal dust emission 
based on the 3D distribution of relativistic electrons and dust grains on the Galaxy, 
and on the 3D pattern of the Galactic magnetic field. The paper is structured as follows: 
Sect.[2]describes the 408 MHz all-sky continnum survey (|Haslam et al 1982|) . the five-year 
WMAP data set ( |Page et al 2007[ ) and the Archeops data (|Ponthieu et al 2005|) used in 
the analysis. In Sect(3]we describe in detail the models. Section 2] describes the statistical 
comparison of data and models. In Sect. [S] we discuss the impact of polarized foreground 
emission on the measurement of the polarized CMB emission with the Planck satellite. 
Finally, conclusions are presented in Section [B] 

2. Observational data 

2.1. Diffuse Galactic synchrotron emission 

The Galactic diffuse synchrotron emission is important at radio and microwave observa- 
tion frequencies. Although its SED is not accurately known it can be well represented by 
a power law in antena temperature l v oc with the synchrotron spectral index ranging 
from —2.7 to —3.3 (K ogut et al 2007[ IGold et al 2009[) . Radio frequency surveys as the 
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Fig. 1. Intensity maps at 408 MHz in K#j units for the Haslam data (left) and built 
with the mode of synchrotron emission with MLS magnetic field for the best fit model 
parameters (right). 



Leider 408 MHz and 1.4 GHz surveys ( jBrouw fc Spoelstra 19761 IWolleben et al 2006|) . 
the Parkes survey at 2.4 GHz (jDuncan 1999ft . and the MGLS survey (Medium Galactic 



Latitude Survey) at 1.4 GHz (Uyaniker et al 1999) are generally used to trace the Galactic 
diffuse synchrotron emission in intensity. For polarization Faraday rotation makes things 
more complex as for frequencies lower than 10 GHz strong depolarization is expected 
(jBurn et al 19661 ISun et al 20081 IJaffe et al 20091 IJansson et al 2009j) so that the best 
Galactic diffuse synchrotron tracers are at high frequency like the WMAP survey at 
23 GHz ( |Page et al 2007] . 



2.1.1. 408 MHz all-sky continuum survey 

In the following we use the 408 MHz all-sky continuum survey (jHaslam et al 1982|) . 
after subtraction of the free-free emission at this frequency, as a tracer of the Galactic 
synchrotron emisison in temperature. We use the HEALpix (IGorski et al 2005)) format 
map available on the LAMBDA website 0. The calibration scale of this survey is claimed 
to be accurate to better than 10 % and the average zero level has an uncertainty of ±3 K as 
explained in lHaslam et al 19821 To substract the free- free emission at 408 MHz we use the 
five-year public WMAP ( Wilkinson Microwave Anisotropics Probe) free-free foreground 
map at 23 GHz generated from the maximum entropy method (MEM) described in 
Hinsha w et al 20071 We have found that the free- free correction has no impact on the final 
results presented on this paper. We start from the full-sky HEALpix maps at N S id e = 512 
(pixel size of 6.9 arcmin) and downgrade them to N S id e = 32 (pixel size of 27.5 arcmin). 
We then substract from the Haslam data the free-free component extrapolated from the 
K-band assuming a power-law dependence of v~ 2A as in IDickinson et al 20031 The left 
panncl of Figure Q] shows the free- free corrected 408 MHz all-sky survey where we clearly 
observe the Galactic plane and the North Celestial Spur at high Galactic latitude. 
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Fig. 2. Form top to bottom: Maps in intensity, I, and polarization Q and U at 23 GHz for 
the WMAP 5-year data (left) and the model of synchrotron emission with MLS magnetic 
field for the best fit model parameters (right) and at 353 GHz for the Archeops data (left) 
and the model of thermal dust emission with MLS magnetic field for the best fit model 
parameters (right). The 353 GHz maps are rotated by 180° for better visualization. All 
the maps are in Krj units. 
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2.1.2. Five-year WMAP polarized data at 23 GHz 



7 



To trace the polarized synchrotron emission we use the all-sky five-year WMAP Q and U 
maps at 23 GHz available on the LAMBDA website in the HEALPix pixelisation scheme 
(jGorski et a l 2005). In particular the coadded maps at N S id e = 512 for the K-band 
( Page et al 2003 IGold et al 2009p . These maps are then downgraded to N si d e = 32 to 
increase the signal-to-noise ratio as we are only interested on very large angular scales and 
the analysis will be performed on Galactic latitude profiles. The pixel window function 
will be accounted for in the following. We assume anisotropic white noise on the maps 
and compute the variance per pixel using the variance per hit provided on the LAMBDA 
website and the maps of number of hits. The large angular scale correlations on the noise 
are neglected. We think this has no impact on the final results as very similar results have 
been obtained using a pixel-to-pixel analysis at N S id e = 16 using the noise correlation 
matrix. The second and third plots on the left column of Figure [5] show the 23 GH Q 
and U maps. We can clearly observe the Galactic plane but also large-scale high Galactic 
structures. 



2.2. Thermal dust 



The dust thermal emission in temperature is well traced by the IRAS (Schlegel et al 1998) 



all-sky observations in the infrared, the FIRAS-COBE (Boulanger et al 1996) all- 
sky observations in the radio and millimeter domain and the ARCHEOPS 
(jMacfas-Perez c t al 2007, Benoit et al 2004]) data in the millimeter domain over 30 % 
of the sky . In polarization early observations by IHiltner 19491 IHall 19491 and later 
by IHeiles 20001 have proved that the starlight emission in the optical domain was polar- 
ized and therefore we will expect the thermal dust emission at millimetric wavelength 
also to be polarized. This was confirmed by the Archeops observations at 353 GHz 
(jPonthieu c t al 2005} that yielded a polarization fraction of about 10 % in the Galactic 
plane. 

Here we use the Archeops 353 GHz Q and U maps as tracers of the polarized thermal 
dust emission. As shown in the fifth and sixth plots of the left column of Figure [5] they 
cover about 30 % of the sky with 13 arcmin resolution. In contrast to the WMAP data at 
23 GHz, most the signal in the Q and U 353 GHz maps is concentrated on the Galactic 
plane. These maps are then downgraded to N S id e = 32 to increase the signal-to-noise 
ratio. As for the WMAP 23 GHz data the noise assumed to be anisotropic white noise 
on the maps and compute the variance per pixel using the variance provided by the 
Archeops collaboration (jMacias-Perez et al 2007[) . 



1 http://lambda.gsfc.nasa.gov/ 
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Fig. 3. Schematic view of the direction of polarization of the Galactic synchrotron and 
dust thermal emission as a function of the Galactic magnetic field direction. 

3. 3D modeling of the galaxy 

We present in this section a realistic model of the diffuse synchrotron and dust polarized 
emission using a 3D model of the Galactic magnetic held and of the matter density in 
the Galaxy. We will consider the distribution of relativistic cosmic-ray electrons (CREs), 
?icre, for the synchrotron emission and the distribution of dust grains, n^nst, for the 
dust thermal emission. The total polarized foreground emission observed at a given 
position on the sky n and at a frequency v can be computed by integrating along the 
line of sight as follows. 



Synchrotron 



For the synchrotron emission (Ribicki & Lightman 1979) and defining 



dJ, 



sync sync 



e svnc (v) n CRE (n,z) 



(5,(n,z) 2 + 5 t (n,z) 2 ) (s+1)/4 dz 



we obtain 
It ync (n) -■ 
Qt vnc (n) -- 
m vnc (n) = 



(1) 



I 

yd/^" c cos(2 7 (n,z)) p sync , 
J d/^" c sin(2 7 (n,z)) p sync . 



(2) 
(3) 
(4) 



where I, Q and U are the Stokes parameters and e sync (v) is an emissivity term. 7 is 
the polarization angle and B n , Bi and B t are the magnetic field components along, 
longitudinal and transverse to the line of sight, n, and z is a ID coordinate along the 
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line-of-sight. s is the exponent of the power-law representing the energy distribution of 
relativistic electrons in the Galaxy. The polarization fraction, p s f™ c i is related to s, as 
follows 

P sync = ( 5 ) 
s + 7/3 

In the following we will assume a constant value of 3 for s so that the synchrotron 
emission will be proportional to the square of the perpendicular component of the 



Galactic magnetic field to the line of sight and p sync = 0.75 (Ribicki fc Lightman 19791 



Locally, the direction of polarization will be orthogonal to the magnetic field lines 
and to the line of sight. Then, the polarization angle 7 is given by 

7(n, s) = - arctan . (6) 



2 \Bf(n,z)-B*(n,z), 
Thermal dust 

For thermal dust emission and defining 

d/^ st (n) = e dust (v) n dust (n,z)dz (7) 

we can write 

I* ust (n) = J Alt us \ (8) 

Qt ust (n) = J d^ tls y ust cos(2 7 (n,z)) 

■/ g (n,z)/ma(n,z), (9) 
U dust (n) = J dl* ust p dust sin(27(n, z)) 

•/ g (n,z)/ ma (n,z). (10) 

where e dust is the dust emissivity, p dust is the polarization fraction, 7 is the polarization 
angle, and, / g and / ma arc polarization suppression factors (see below). Hereafter, the 
polarization fraction p dust is set to 10% following IPonthieu et al 20051 



As discussed before, dust grains in the ISM arc oblate and will align with their large 
axis (see Figure [3]) perpendicularly to the magnetic field lines (jDavis fc Greenstein 19511 
ILazarian 19951 lLazarian 19971 Lazaria n 200"9|) . Therefore, the polarization direction, for 
thermal dust emission will be perpendicular both to the magnetic field lines and the 
line-of-sight as was already the case for the synchrotron emission. Then, the polarization 
angle 7 will be the same for the synchrotron and dust emission. However, as the dust 
grains rotate with their spin axis parallel to the magnetic field we also need to account 
for a geometrical supression factor. For example, if the magnetic field direction is parallel 
to the line-of-sight we expect the dust polarized emission to be fully supressed. The 
suppression factor can be expressed as / g = sin 2 (a) where a is the angle between the 
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magnetic field direction and the linc-of-sight. By construction, we observe that 7 and a 
are the same angle. The process of aligncmcnt of the dust grains with the magnetic held 
is very complex (jMathis 19861 IGoodman fc Whittet 19951 ILazarian 19951 ILazarian 19971 
Lazaria n 2009[) and its accurate representation is out of the scope of this paper. Then, 
to account for missalignements between the dust grains and the magnetic held lines we 
dehne an empirical factor f ma . The form of this factor is unknown but we have empirically 
observed that the geometrical supression seems to be more important than expected for 
the Archeops data. Therefore, we have taken f ma to be oc sin(a). We have observed that 
the results presented on this paper are robust with respect to this parameter. 



3.1. Matter density model 

In galactocentric cylindrical coordinates (r, z, <f) the relativistic electron density distri- 
bution can be written as (see |Drimmel fe; Spergel 2001 1 



«.CReO,2) = n 0>e r-—r, r 7-—, r, (11) 

cosh(z/kpc) cosh(z/ncRE,/i) 
where ncRE,/i defines the width of the distribution vertically and is setted to 1 kpc 
in the following. ncRE,r defines the distribution radially and it is a free parameter 
of the model. Notice that we expect these two paramters to strongly correlated and 
that is why we have decided to fix one of them as was the case on previous analyses 
(jSun et al 20081 IJaffe et al 2009p . 

The density distribution of dust grains in the Galaxy is poorly known and that is why 
we have chosen to describe it in the same way that the relativistic electrons one 

n d (r,z) = n 0td -, (12) 

cosh (z/nd,h) 

where n d ^ r and nd,h are the radial and vertical widths of the distribution. In the following 
we set them to 3 and 0.1 kpc respectively. 



3.2. Galactic magnetic field model 

According to observations, many spiral galaxies in our vicinity and for var- 
ious redshifts, present a large scale magnetic field with intensity of the 
order of few fiG and direction spatially correlated with the spiral arms 
(|Sofue et al 19861 IBeck et al 19961 IWielebinski 20051) . For our Galaxy, the mag- 
netic field direction also seems to follow the spiral arms but with a complex 
spatial distribution (IWielebinski 20051 Man et al 20061 IBeck 2006j) .Indeed, there 
are hints for local inversion of the field direction and radial dependency of 
the intensity (B an et al 20061 IBeck 200"Tj) . Pulsar Faraday rotation measurements 
Pan et al 20041 IHan et al 20061 ISofue et al 19861 IBrown et al 2007j) have been used to 
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fit the Galactic large-scale magnetic field with various models including a axi-symmetric 
or bi-symmetric form, a field that reverses in the inter-arm regions, etc. Pulsar Faraday 
rotation measurements also indicate the presence of a turbulent component of the 
magnetic field ([Han et al 2004p . 



3.2.1. Large-scale magnetic field. 

In the following we consider a Modified Logarithmic Spiral (MLS) model of the large- 
scale magnetic field based on the WMAP team model presented in Page et al (2007). 
It assumes a logarithmic spiral to mimic the shape of the spiral arms (jSofue et al 1986|) 
to which we have added a vertical component. In galactocentric cylindrical coordinates 
(r, z, $) it reads 



B(r) = S reg (r)[cos(0 + /3)ln f^-J sin(p)cos(x) • u r 
- cos((f> + 13) In cos (p) cos (x) ' u <t> 

+sin(x) ■ u z ], (13) 

where p is the pitch angle and (3 = 1/ tan(p). r$ is the radial scale and x{ r ) = Xo{ r ){z y ' Zq) 
is the vertical scale. Following |Tay lor fc Cordes 1993 we restrict our model to the range 



3 < r < 20 kpc. The lower limit is setted to avoid the center of the Galaxy for which the 
physics is poorly constrained and the model diverges. The intensity of the regular field 
is fixed using pulsar Faraday rotation measurememts bv lHan et al 20061 

Breg{r) = B e" 1 ^ (14) 

where the large-scale field intensity at the Sun position is Bq = 2.1 ± 0.3/iG and 
Rb = 8.5 ± 4.7fcpc. The distance between the Sun and the Galactic center, R@ is set to 
8 kpc (jEisenhauer et al 20031 fReid fc Brunthaler 2005]) . 

Wc also study the spiral model of IStanev 19971 ISun et al 20081 hereafter ASS. In 
cylindrical coordinates it is given by 

B° = D 1 (r,$,z)D 2 (r,<?>,z)sin(p) (15) 
B§ = -D 1 (r,^,z)D 2 (r,^,z)cos(p) (16) 
B? = (17) 

where D\ accounts for the spatial variations of the field and D 2 for asymmetries or 
reversals in the direction. The pitch angle is defined as for the MLS model described 
above. Di(r, z) is given by 
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Di(r,z) 




[ ) r>R c 
r<R c 



(18) 



where Rq is the distance of the Sun to the center of the Galaxy and it is set to 8 kpc as 
before. R c is a critical radius and it is set to 5 kpc following the ASS+RING model in 
ISun et al 20081 In the same way Rq is fixed to 10 kpc, Bo to 6 /J.G and B c to 2 /j,G. The 
field reversals are specified by 



D 2 (r) 



-1 r>7.5 

-1 6<r<7.5 

-1 5 < r < 6 

-1 r < 5 



(19) 



where r is defined in units of kpes. 

It is important to notice that the synchrotron and thermal dust polarized emission de- 
pends only on the orientation and not on the sign of the magnetic field and therefore, 
they are not sensible to field reverses in the model. 



3.2.2. Turbulent component 

In addition to the large-scale Galactic magnetic field, Faraday rotation measurements on 
pulsars in our vicinity have revealed a turbulent component for scales smaller than few 
hundred pc ( |Lyne fc Smith 1989| ) . Moreover this turbulent field seems to be present at 
large scales (|Han et al 2004[) . The amplitude of this turbulent component is estimated to 
bc of the same order of magnitude as the amplitude of the regular one (|Han et al 2006[) . 
The magnetic energy i?e(fc) associated with the turbulent component is well described 
by a Kolmogorov spectrum (|Han et al 20041 Irian et al 2006[) 

E B (k) = C (£f (20) 

where a = -0.37 and C = (6.8 ± 0.3) ■ 10~ 13 erg cm" 3 kpc. 

3.2.3. Final model 

Finally the total magnetic Galactic field £? tot (r) can be written as 

B tot (v) = B reg (v) + B turb (r) (21) 

where B reg (r) is the regular component, either MLS or ASS, and B tur b(r) is the 
turbulent one. We define A tur b as the relative intensity of the turbulent component 
with respect to the regular one and it is a free parameter of the model. The turbulent 
component is computed from a 3D random realization of the Kolmogorov spectrum 
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In this paper we do not consider the halo component presented by ISun et aT~2 008 , 
IJansson et al 2009f as relativistic electrons and dust grains are not expected to be con- 
centrated on the halo. 



3.3. Emissivity model in polarization 

As discussed in the previous section, the polarized emission in the 23 GHz WMAP data 
shows complex structures both on the Galactic plane and in local high Galactic latitude 
structures such as the North Celestial Spur ( |Page et al 2007[ ) . An accurate representation 
of this complexity can not be achieved using a simplified model like the one presented 
here. A similar degree of complexity is observed in the 353 GHz polarization maps al- 
though the morphology of the structures is rather different. To account for this complexity 
the Q and U estimate for synchrotron and thermal dust models are corrected using in- 
tensity templates of these components extrapolated to the observation frequencies (23 
and 353 GHz) using constant spectral indeces. 
For the synchrotron emission we have 



Qs - /Has [ QAm ) jsync , (22) 

/ v \P> U* ync 

U s = ^Has L ,„„ "sync ) (23) 



where -?Has is the reference map in intensity constructed from the 408 MHz all-sky sur- 
vey (section I2.1.1|) after subtraction of the free-free emission and v is the frequency of 
observation. Notice that we do not use the synchrotron MEM intensity map at 23 GHz 
(jHinshaw et al 2007[) as a synchrotron template to avoid any possible spinning dust con- 
tamination. The spectral index (3 S used to extrapolate maps at various frequencies is 
assumed to be spatially constant on the sky and is a free parameter of the model. 
For the thermal dust emission we write 

Qdust 

Q d = Is f d -nks7> ( 24 ) 

jjdust 

U d = ^fdjtlT, (25) 

where I s fd is the reference map in intensity at 353 GHz generated using model 8 
of Schlegel, Finkbeiner and Davis (jFinkbcincr et al 1999|) which was obtained from the 
IRAS [Schlegel et al 1998] and COBE/DIRBE data ( jBoulanger et al 1996[ ). 

We compute the above I,Q and U maps for synchrotron and thermal dust with the 
help of a modified version of the Hammurabi code ([Waclkcns ct al 2009]). Each map is 
generated by integrating in 100 steps along each line-of-sight defined by the HEALPix 
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Table 1. Latitude and longitude bands for the Galactic profiles used in the analysis. 



Latitude interval (deg) 


[0,30] 


[30, 90] 


[90, 120] 


[120, 180] 


[180, 270] 


[270, 330] 


[330, 360] 


Longitude interval (deg) 


[-90, -50] 


[-50, -20] 


[-20, -5] 


[-5,5] 


[5,50] 


[50, 70] 


[70, 90] 



Table 2. Parameters of the 3D Galactic model. 



Parameter 


Range 


Binning 


P (deg) 


[-80.0,80.0] 


10.0 


-^turb 


[0, 2.5] * B reg 


0.25 


flCRE.r (kpc) 


[0.0, 10.0] 


1 


Pa 


[-4.3, -2.4] 


0.1 



•Nsidc = 128 pixel centres. The integration continues out to 25 kpc from the observer 
situated 8.5 kpc from the Galactic centre. 

4. Galactic-profiles comparison 

4.1. Galactic-profiles description 

In order to compare the models of Galactic polarized emission to the available data 
we compute Galactic longitude and latitude profiles for the models and for the data in 
temperature and polarization using the set of latitude and longitude bands defined in 
Table Q] For the Galactic longitude and latitude profiles both, we use bins of longitude of 
2.5°. In the following discussions we will only consider Galactic latitude profiles because 
equivalent results are obtained with the longitudinal profiles. 

We compute error bars including intrinsic instrumental uncertainties and the extra 
variance induced by the presence of a turbulent component. The latter is estimated from 
the RMS within each of the latitude bins following [Jansson et al 20091 For the 408 MHz 
all-sky survey we account for intrinsic uncertainties due to the 10 % calibration errors 
described in Section [5] For the WMAP 23 GHz data we have computed 600 realizations 
of Gaussian noise maps from the number of hits per pixel and the sensitivity per hit 
given on the WMAP LAMBDA web site. We have computed Galactic latitude profiles in 
polarization for these simulated maps and estimated intrinsic errors from the standard 
deviation within each latitude bin. For the Archeops data we use the noise simulations 
discussed in IMarias-Perez et al 20071 and proceed as for the WMAP data. 
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Fig. 4. Galactic profiles in temperature at 408 MHz built using the Haslam data (black) 
and the model of synchrotron emission with MLS for various values of the pitch angle p 
(form green to red). 



Galactic latitude profiles are computed for a grid of models obtained by varying 
the pitch angle, p, the turbulent component amplitude, A tur b, the radial scale for the 
synchrotron emission, ncRE.r and the synchrotron spectral index, j3 s . The range and 
binning step cosidered for each of these parameters is given in Table [2 All the other 
parameters of the model of the Galactic magnetic field and matter density are fixed to 
values proposed in Section EH Notice that to be able to compare the dust models to the 
Archeops 353 GHz data, the simulated maps are multiplied by a mask accounting for 
the Archeops incomplete sky coverage. 
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Fig. 5. Galactic profiles in polarization Q and U at 23 GHz built with the five-year 
WMAP data (black) and the model of synchrotron emission with MLS magnetic field for 
various values of the pitch angle p (form green to red). 



Figure E] shows in black Galactic latitude profiles in temperature for the 408 MHz all- 
sky survey with error bars computed as discussed above. In colors, we show for comparison 
the expected synchrotron diffuse Galactic emission from the MLS Galactic magnetic field 
model for various values of the pitch angle p. On Figure [5] we present the polarization 
Galactic latitude profiles for the WMAP 23 GHz data (black) and the expected polar- 
ized diffuse synchrotron emission for the previous MLS models (color). Finally, Figure |5] 
shows the polarization Galactic latitude profiles for the 353 GHz Archeops data (black) 
compared to the same MLS models (color). From these figures we can see that the cur- 
rent available data do help discriminating between the different models and therefore a 
likelihood analysis is justified. 
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Fig. 6. Galactic profiles in temperature and polarization Q and U at 353 GHz with the 
ARCHEOPS data (black) and for various values of the pitch angle p for the model in 
polarization of thermal dust emission with MLS magnetic field (from green to red). 

4.2. Likelihood analysis 

The data and model Galactic latitude profiles are compared using a likelihood analysis 
where the total likelihood function is obtained from 

dot = nLiAi (26) 

where for each of the 3 data sets described above the log-likelihood function is given by 

Ni on -l N lat -l (n d _ nrd \2 

-iog^ = E E E , * ( 27 ) 

where i represents the polarization state meaning intensity only for the 408 MHz 
all-sky survey, and, Q and U polarization for the 23 GHz WMAP and 353 GHz 
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Table 3. Best-fit parameters for the MLS and ASS models of the Galactic magnetic 
field. 



Data 


Magnetic field model 


p(deg) 




HCRE.r 




Xmin 


408 MHz 


MLS 


-20.ot|° : ° 


< 1.00 (95.4 % CL) 


4±^ 6 





3.58 


ASS 


-io.otr 


< 1.0 (95.4 % CL) 


c + 15 
°-3 





4.65 


WMAP 23 GHz 


MLS 


-30 0+ 40 


< 1.25 (95.4 % CL) 


< 20 (95.4 % CL) 


_o 4 +o.i 


5.72 


ASS 


-40 0+ 60 - 

^U.U_3q.O 


< 1.5 (95.4 % CL) 


31^(95.4% CL) 


_3 4+ - 1 

°-^-0.8 


7.62 


Archeops 353 GHz 


MLS 


-20+ 80 


< 2.25(95.4%CL) 








1.98 


ASS 


60.0±lo 


25+ 2 








1.72 


All 


MLS 


-20i 8 5 ° 


< 2.25(95.4%CL) 








1.98 


ASS 


60.0+4°, 


25+ 2 








1.72 



Archeops data, j and k represents the longitude bands and latitude bins respectively. 
Dfj k and Mfj k corresponds to the data set d and model for the i polarization state, j 
longitude band and k latitude bin, respectively, afj k is the error bar associated to . k . 



Table [3] presents the best-fit parameters for the three-data sets described above and 
also for the combination of three of them (labeled All in the table). Results are pre- 
sented both for the MLS and ASS models of the Galactic magnetic field. The best-fit 
value for the pitch angle, p, are in agreement within 1-cr error bars for the three data 
sets but for the Archeops data in the case of an ASS model. Notice that our results 
are compatible with the pitch angle values presented in ISun et al 20 08 , Page et al 2007| 
Mivillc-Dcschcnes et al 2008] The relative amplitude of the turbulent component, A tur b, 
is poorly constrained and the data does not seem to favour a strong turbulent com- 
ponent either in the case of MLS or ASS models. However, our results are compatible 
with the ones presented in ISun et al 20081 IMiville-Deschenes et al~20 08. Han et al 20041 
Ha n et al 20061 at the 2-cr level. The electronic density radial scale, noRE.rj is poorly con- 
strained by the data both for MLS and ASS models although our results are compatible 
with those of ISun et al 20081 We also tested the possibility of a local contribution to 
the electronic density as proposed by ISun et al 20081 We found that adding this local 
component does not improve either the fit or the constraint on the radial scale. The 
best-fit value for the spectral index of the synchrotron emission seems to be significantly 
lower than the one in ISun et al 200l |Page et al 2007] This may due to differences in the 



intensity template. Notice that we rescale the polarization intensity using the 408 MHz 
all-sky survey to obtain a more realistic model. 
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Fig. 7. From left to right and from top to bottom : power spectra 
CP \Cf E \Cf B \Cj E \Cj B \Cf B at 23 GHz built with the WMAP 5-year data (black) 
and the model of synchrotron emission with MLS magnetic field for the best fit model 
parameters, excluding the Galactic region defined by |b| < 5°. 



4.3. Temperature and polarization angular power spectra 

Using the best-fit parameters of the MLS model, p = —30.0°, A tur b = 0.0, ncKE,r = 1 
and (3 S = —3.4, we have constructed simulated maps of the sky at 408 MHz and 23 and 
353 GHz. These maps are shown on right-hand side of Figures 12.1.11 and [2] Although 
for the Galactic profiles the fit can be considered relatively good, the fake temperature 
map at 408 MHz looks very different from the 408 MHz all-sky survey map (left side 
of the plot), in particular at the North Polar Spur (fWollcbcn et al 2007[) . as no local 
structures were included in the model. This supports a posteriori our correction of the 
polarization synchrotron model using an intensity template as presented in Section [3] 
In Q and U polarization, the 23 GHz fake maps seem to reproduce qualitatively the 
structure observed in the WMAP data (left side of the plot). However in temperature the 
model and the data are very different as we have not account for a variable synchrotron 
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Fig. 8. From left to right and from top to bottom: power spectra 
CP \Cf E \Cf B \Cj E \Cj B \Cf B at 353 GHz computed from ARCHEOPS data (black) 
and the model of thermal dust emission with MLS magnetic field for the best fit model 
parameters (red) for the full-sky. 

spectral nor fo any extra component as discussed in |Page et al 2007| |Kogut et al 2007} 
IMiville-D eschenc s et al 20081 [Macias- Perez et al 20101 Finally, the model of thermal 
dust emission is able to reproduce qualitatively the Archeops data at 353 GHz. 

Figures [7] and [5] show the temperature and polarization angular power spectra 
for the 23 GHz WMAP and 353 GHZ Archeops data compared to the best-fit MLS 
model for synchrotron and dust, respectively. As discussed before, the temperature 
auto power spectrum of the 23 GHz data is very different from the model as no 
extra components in temperature were considered. However in polarization we have 
qualitatively a good agrement. However, we clearly observe that the model does not 
account for all the observed emission. At 353 GHz the agreement between the data 
and the model qualitatively and quantively is good. For polarization most of the data 




Fig. 9. From left to right and from top to bottom: power spectra 
c tt £EE c bb £TE qTB c eb at 10Q GHz the model of Ga i actic polarized emis- 
sion applying a Galactic cut of \b\ < 15° (black) 30° (blue) and 40° (cyan) and 
compared them to the expected CMB ones (red) for the WMAP best-fit ACDM model 
(jKomatsu et al 2008[) to which we added a tensor component assuming an tensor to 
scalar ratio of 0.1.. 



samples at less than 3-cr from the model. In temperature the model is not as accurate 
as in polarization but we notice that the fitting was restricted to polarization data only. 



5. Galactic foreground contamination to the CMB measurements by the 
Planck satellite 

We can use the best-fit model of polarized synchrotron and dust emission to estimate 
the polarized foreground contamination to the CMB at the Planck satellite observation 
frequencies. We are interested in comparing the expected foreground contribution to the 
expected CMB emission. Notice that the aim of this section is not to obtain an accurate 
template of the polarized Galactic foreground emission to be subtracted on the Planck 
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data for CMB purposes. 

For this purpose, we have produced simulated maps of the Galactic polarized fore- 
ground emission using the best-fit model parameters for each of the Planck CMB fre- 
quencies, 70, 100, 143 and 217 GHz. The thermal dust polarized emission have been 
extrapolated using a constant spectral index of 2.0 in antenna temperature. We have 
computed the temperature and polarization power spectra of these maps and compared 
them to the expected CMB ones for the WMAP best-fit ACDM model (available on the 
LAMBDA website) to which we added a tensor component assuming a tensor to scalar 
ratio of 0.1. Note that neither noise nor systematics nor resolution effects are considered. 

Figure O shows these power spectra at 100 GHz. The expected CMB signal is rep- 
resented in red. The polarized diffuse foreground emission for Galactic latitude cuts of 
|6| < 15°, 30° and 40° are shown in solid black, blue and cyan lines, respectively. 1-a 
errors in the model are represented are dashed lines. In temperature, the CMB Cf T 
dominates at all the angular scales considered as could be expected from the WMAP 
and Archeops data. For polarization, the CMB C EE dominates at high I values but we 
observe significant foreground contamination at the lowest i values [t < 20). In the same 
way, the CMB Cj E dominates at 100 GHz but for very low t values. However, the CMB 
C BB is significantly smaller than the foreground contribution at all the angular scales 
considered even for such a large value of the tensor to scalar ratio. The CMB Cf B and 
C EB are expected for most cosmological models to be null and therefore, the foreground 
contribution dominates the signal. 

As the Galactic polarized foreground emission seems to dominate the observed emis- 
sion at the Planck CMB frequencies, special care should be taken when estimating the 
CMB emission using standard template subtraction techniques and component separa- 
tion algorithms. The assesmcnt of the final errors is crucial and we think that models of 
the polarized foreground emission such as those presented in this paper can be of much 
help for this task. 

6. Summary and conclusions 

We have presented in this paper a detailed study of the diffuse synchrotron and thermal 
dust polarized Galactic foreground emission components. We have constructed coherent 
models of these two foregrounds based on a 3D representation of the Galactic magnetic 
field and of the distribution of relativistic electrons and dust grains in the Galaxy. For the 
Galactic magnetic field we have assumed a large-scale regular component plus a turbulent 
one. The relativistic electron distribution and dust grain distribution have been modeled 
as exponentials peaking at the Galactic center. From these analyses we have been able 
to study the main parameters of the models, the magnetic field pitch angle, p, the radial 
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width of the relativistic electron distribution, h er , the relative amplitude of the turbulent 
component, A tur b and spectral index of the synchrotron emission j3 s . We have been able 
to set constraints only on the pitch angle and the synchrotron spectral index. An upper 
limit on the relative amplitude of the turbulent component is obtained although the data 
seems to prefer no turbulence at large angular scales. With the current data we are not 
able to constrain the radial width of the relativistic electron distribution. Notice that our 
constraints are compatible with those in the literature. 

Using the best-fit parameters we have constructed maps in temperature and polar- 
ization for the synchrotron and dust thermal emission at 23 and 353 GHz and compare 
them to the WMAP and ARCHEOPS data at the same frequencies. We find good agree- 
ment between the data and the model. However, when comparing the temperature and 
polarization power spectra for the data and model maps, we observe that synchrotron 
emission model is not realistic enough. For dust the model seems to reproduce better the 
data but it is important to realize that the errors on the Archeops data are much larger. 

From this, we can conclude that the models presented in this paper can not be used 
for direct subtraction of polarized foregrounds for CMB purposes. However they can be 
of great help for estimating the impact of the polarized Galactic foreground emission on 
the reconstruction of the CMB polarized power spectra. Indeed, we have extrapolated 
the expected polarized Galactic foreground emission to the Planck CMB frequencies, 
70, 100, 143 and 217 GHz and found they dominate the emission at low £ values where 
the signature in the polarized CMB power spectra of important physical processes like 
rcionization is expected. Furthermore, the Galactic polarized foreground emission seems 
to dominate the B modes for which we expect an unique signature from primordial 
gravitational waves. Because of this, we propose the use of models like the ones presented 
in this paper to asses the errors in the reconstruction of the CMB emission when using 
template subtraction techniques or component separation algorithms. 
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